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Abstract 

The h c i}P\) state of charmonium has been observed in the reaction ip(2S) — > ir h c — > (77) (777c) 
using 3.08 million tp(2S) decays recorded in the CLEO detector. Data have been analyzed both 
for the inclusive reaction, where the decay products of the n c are not identified, and for exclusive 
reactions, in which rj c decays are reconstructed in seven hadronic decay channels. We find M(h c ) = 
3524.4 ± 0.6 ± 0.4 MeV which corresponds to a hyperfine splitting AM hf (lP) = (M( 3 Pj)) - 
M( 1 P 1 ) = +1.0 ± 0.6 ± 0.4 MeV, and B(ip(2S) -» n°h c ) x B(h c -► 7%) = (4.0 ± 0.8 ± 0.7) x 10~ 4 . 



2 



Over the past thirty years charmonium spectroscopy has provided valuable insight into 
the quark-antiquark interaction of quantum chromodynamics (QCD). QCD-based potential 
models have been quite successful in predicting masses, widths, and dominant decays of 
several charmonium states. The central potential in most of these calculations is assumed 
to be composed of a vector Coulombic potential (~l/r) and a scalar confining potential 
(~r). Under these assumptions, the spin-spin interaction in the lowest order is finite only for 
L = states. It leads to the hyperfine splittings AM h f{nS) = M(n 3 Si)— M(n 1 S , ) between 
spin-triplet and spin-singlet 5-wave states of charmonium, which have been measured as 
AM hf (lS) = M(J/ip)-M(r) c ) = 115 ±2 MeV |lj, AM hf (2S) = M(if>(2S))-M(rf c ) = 
43+3 MeV 00. It also leads to the prediction that the hyperfine splitting AM hf ((M( 3 Pj))- 
M( l Pi)) for P-wave states should be zero. Higher-order corrections are expected to provide 
no more than a few MeV deviation from this result 0,0,0. Lattice QCD calculations 
predict AM^ (IP) = +1.5 to +3.7 MeV, but with uncertainties at the few-MeV level. Larger 
values of AM^(IP) could result if the confinement potential had a vector component or if 
coupled channel effects were important. In order to discriminate between these possibilities, 
it is necessary to identify the h c Q-P\) state and to measure its mass to (9(1 MeV) as the 
mass of the 3 Pj centroid is very well known, (M( 3 Pj)) = 3525.36+0.06 MeV 0, 

In this Letter we report the successful identification of h c in the isospin-violating reaction 

e+e- -> if>(2S) -> n°h c , h e -> jrj c , ^ -> 77- (1) 

Two methods are used: one in which the r] c decays are reconstructed (exclusive), which 
has an advantage in signal purity, and the other in which the r\ c is measured inclusively, 
which has larger signal yield. Together these approaches provide a result of unambiguous 
significance, and allow a precise determination of the mass of h c and the branching fraction 
product B^Bh, where B^ = B{jp{2S) — > n h c ) and B h = B{h c — > JT] C ). Theoretical estimates 
of the product B^Bh vary by nearly two orders of magnitude, (0.5 — 40) x 10~ 4 00. 

The Crystal Ball Collaboration at SLAC searched for h c using the reaction of Eq. (JTJ 
but were only able to set a 95% confidence upper limit B^Bh < 16 x 10~ 4 in the mass 
range M(h c ) = (3515 - 3535) MeV 0. The FNAL E760 Collaboration searched for h c 
in the reaction pp — ► h c — > n^J/ip, J/ip — * e + e~, and reported a statistically significant 
enhancement with M(h c ) = 3526.2+0.15+0.2 MeV, T(h c ) < 1.1 MeV 0. The measurement 
was repeated twice by the successor experiment E835 with ~2 x and ~3 x larger luminosity, 
but no confirming signal for h c was observed in h c — > ix°J/ip decay 0. 

A data sample consisting of 3.08 x 10 6 '0 (25') decays was obtained with the CLEO III 
and CLEO-c detector configurations [13, H3, E3, 03 at the Cornell Electron Storage Ring. 
The CLEO III detector features a solid angle coverage for charged and neutral particles of 
93%. The charged particle tracking system achieves a momentum resolution of ~0.6% at 
1 GeV, and the calorimeter photon energy resolution is 2.2% for P 7 = 1 GeV and 5% at 
100 MeV. Two particle identification systems, one based on energy-loss (dE/dx) in the drift 
chamber and the other a ring imaging Cherenkov (RICH) detector, are used to distinguish 
pions from kaons. 

Half of the ip(2S) data were accumulated with a newer detector configuration, CLEO- 
c |lj|, in which the silicon strip vertex detector was replaced with an all-stereo six- layer 
wire chamber. The two detector configurations also correspond to different accelerator 
lattices. Studies of Monte Carlo simulations and the data reveal no significant differences 
in the capabilities of the two detector configurations; therefore the CLEO III and CLEO-c 
datasets are analyzed together. 
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FIG. 1: Exclusive analysis: (a) Scatter plot of the reconstructed rj c candidate mass vs. the recoil 
mass against ir° for data. The horizontal band in the vicinity of M(J/ip) and the diagonal band 
at larger rj c candidate mass correspond to tp(2S) — > Tr 7r°J/ip and tp(2S) — ► 7Xc0> respectively. 
The dashed lines denote the region M(rj c ) = 2982 ± 50 MeV. Data events (open histograms) 
and Monte Carlo background estimate (shaded histograms) of (b) reconstructed r\ c candidate mass 
projection for M(ir° recoil) = 3524 ±8 MeV and (c) recoil h c candidate mass spectrum for M(r] c ) = 
2982 ±50 MeV. The peaks in b) near M=3.1 GeV and 3.25 GeV correspond to ^(25) -> Tr°7r°J/^ 
and ip(2S) — > 7X c o, respectively, (d) Reconstructed rj c candidate mass for data in the direct decay 
ip(2S) —> -)rj c . The peak near M = 3.4 GeV is from the direct decay ifi(2S) — > 7Xc0- 



The inclusive and exclusive analyses share a common initial sample of events and nu- 
merous selection criteria. Details of the analyses are provided in a companion paper ^ . 



Event selection for both analyses require at least three electromagnetic showers and two 
charged tracks, each selected with standard CLEO criteria. For showers, £7 7 > 30 MeV is 
required. Candidates for 77 decays of tt° or rj mesons satisfy the requirement that M(77) be 
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FIG. 2: Inclusive analysis: Recoil mass against 7r° for (a) Monte Carlo sample for 39.1 million 
ip{2S) and (b) data for 3 million tp(2S). See text for details. 



within 3 standard deviations (a) of the known tt° or rj mass, respectively. These candidates 
are kinematically fit, constraining M^jj) to the appropriate mass to improve 7r°/r) energy 
resolution. Charged tracks are required to have well-measured momenta and to satisfy cri- 
teria based on the track fit quality. They must also be consistent with originating from the 
interaction point in three dimensions. 

Both techniques identify h c as an enhancement in the spectrum of neutral pions from the 
reaction ip(2S) —> ix°h c 15;]. For this purpose, it is useful to remove neutral pions originating 
from any other reaction. It is easy to remove most of the ir° arising from ip(2S) — > 7t + tt~ J/ip, 
with J/ip —>■ 7T° + hadrons and J/ip, with J/ip — > any. Since the recoil spectra against 
M(tt + tt~) (both analyses) and M(7r°7r°) (inclusive only) show prominent peaks for J/ip] 
these events are removed by appropriate selection around M(J/ip). 

In the exclusive analysis, 7] c are reconstructed in seven channels: K s K ± it T , K^K^n^ , 



K + K~7t + tt~ , 7r + 7T~7r + 7r~, K + K~it°, TT + 7r~7](^ 77), and 7T + 7T~ri(^ 7r + 7T _ 7r°). The sum of 
the branching fractions is (9.7 ± 2.7)% 0. The decay chain in Eq. as well as these r\ c 
decays are identified from reconstructed charged particles, n° and r\ mesons. For r] decays to 
7r + 7T~7r°, the three-pion invariant mass is required to be within 20 MeV of the nominal f] mass. 
The Kg candidates are selected from pairs of oppositely charged and vertex-constrained 
tracks with invariant mass within 10 MeV, roughly 4<r, of the Kg mass. A kinematically 
constrained 4C fit is performed for each event. A 1C fit is performed for the i] c — ► K^K ± n T 
decay because the K Q L is not detected. It is required that M(r] c ) = 2980 ± 50 MeV. No 
explicit selection of the energy of the photon from h c — > 7% is required. The final selection 
is on the r] c candidate mass; however, to improve resolution, the h c mass is calculated from 
the 4-momentum of the if)(2S) and the 7r° instead of the invariant mass of its decay products. 

In addition to ip(2S) — > nnJ/ip decays discussed above, a fraction ofip(2S) decays proceed 
through ip{2S) — > 7r°J/^ and ip(2S) JXcj ~^ tt°X. To suppress the 7r° background, each 
signal photon candidate is paired with all other photons in that event. If the invariant mass 
of any pair is within the 7r° mass requirement, the event is removed. 
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Fig. E^a) shows the scatter plot of the r\ c candidate mass versus 7r° recoil mass (sum of 
all channels). Many events are seen in the vicinity of M(J/ip). In the mass band M(j] c ) = 
2980 ±50 MeV an enhancement of events is observed at larger 7r° recoil mass. The projection 
of the events in this band and the Monte Carlo background estimate is shown in Fig. ^c). 
A prominent peak is clearly visible over a very small background. The projection of the 
events in the mass band M(ir Q recoil) = 3524 ± 8 MeV and the Monte Carlo background 
estimate, shown in Fig. ^b), indicate that most of these events arise from r] c decay. The 
7T° recoil mass spectrum, in Fig. H^b), is fit using a double Gaussian shape determined from 
Monte Carlo simulation and an ARGUS function background 0] . The maximum likelihood 
fit yields 17.5 ± 4.5 counts in the peak and M(h c ) = 3523.6 ± 0.9 MeV. 

Several different methods have been utilized to estimate the statistical significance s 
of the signal including the fit to the recoil mass spectrum just described, Poisson 
fluctuations of MC-predicted backgrounds inside the signal window, and a binomial statistics 
calculation using the assumption that the events in the recoil mass distribution are uniformly 
distributed. Using the difference between the likelihood values of the fit with and without 
the signal contribution, we obtain s = 6. lex; similar calculations with different rj c mass ranges 
yield s = 5.5 — 6.6a. The probability that Poisson fluctuations of the background, estimated 
from the generic MC sample, completely account for the observed events in the signal region 
is 1 x 10~ 9 (s = 6.0a). The binomial probability that the number of data events in Fig. ^b) 
and^c) fluctuate to be greater than the number of events in the signal region is 2.2 x 10~ 7 , 
corresponding to s ~5.2<r. 

To test our ability to reconstruct i] c decays and provide normalization for the branching 
fraction measurement, B^Bh, the direct radiative decay ip{2S) — > jr/c is studied. Events 
are reconstructed in the same r) c decay channels as for the h c search, but with much better 
yields. Relative yields among the various channels are similar to previous results and the 
1] C peak shape was verified for each channel. Figures ^b) and E^d) show the reconstructed 
mass spectra for the rj c candidates from h c and direct ip(2S) decay, respectively. The rj c 
mass resolution in the photon recoil mass spectrum is identical for all seven channels. This 
distribution summed over all channels (not shown) is fit using a peak shape which consists 
of a Monte Carlo-derived double Gaussian convolved with a Breit-Wigner function (with 
M(rj c ) = 2979.7 MeV, F(t] c ) = 27 MeV). It yields 220 ± 22 counts. The efficiency-corrected 
ratio of h c decays to direct decays, which corresponds to ^i^/i^jwhere £>© = B(ip{2S) — > 
7^ c ), is determined to be 0.178±0.049. The CLEO [13 and PDG pj values are combined to 
obtain B v = (0.296±0.046)%. Multiplying these two results yields B^B h = (5.3±1.5)xl0" 4 
from the exclusive analysis. 

In the inclusive analysis, we explore two methods to enhance the selection of neutral 
pions which are part of the chain ip{2S) — > 7r°h c — > 7T ^i] c . One way is to specify that there 
be only one photon in the event with energy for the transition h c — > / ~fi] c , fa 503 MeV 
(corresponding to M(h c ) pa 3526 MeV). Another way is to specify that the mass recoiling 
against the photon and 7r° for the event should be near the mass of i] c . Both approaches are 
investigated, leading to consistent results, as detailed in Ref. 

A combined sample of generic ip{2S) decay and signal Monte Carlo events is used to 
optimize the criteria for the final event selection. The resulting selection criteria determined 
were E 1 = 503 ±35 MeV for hard photon acceptance in one approach and M(rj c ) = 2980 ±35 
MeV in the other. As a result of the Monte Carlo studies, a number of selection criteria, in 
which the two approaches occasionally differ, are made. These include requiring only one 
7T° in the signal region, removing hard photons that reconstruct 7] mesons with any other 
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FIG. 3: Inclusive analysis: Efficiency-corrected fit yields versus | cos#| for data with E 1 = 503 ± 
35 MeV, (a) for h c yield, the curve corresponds to the best fit oc (1 + cos 2 6) and (b) for the nearly 
isotropic background yield. 



TABLE I: Results for the inclusive and exclusive analyses for the reaction ip(2S) — > ir°h c — ► 7r°777 c . 
First errors are statistical, and the second errors are systematic, as described in the text and 
Table HH 



Inclusive Exclusive 



Counts 150±40 17.5±4.5 

Significance ~3.8er 6.1a 

M{h c ) (MeV) 3524.9±0.7±0.4 3523.6±0.9±0.5 

B^B h (10~ 4 ) 3.5±1.0±0.7 5.3±1.5±1.0 



photon, accepting photons in the calorimeter endcaps, removing photons from the cascade 
reaction ip(2S) — > 7Xj — > llJ/^i an d the choice of the background shape. 

The recoil spectrum for the total Monte Carlo sample of 39.1 million ip(2S) (13 times the 
size of the data sample), obtained in the _E 7 -seIection approach with its optimized selection 
criteria, is shown in Fig. Efa)- A product branching fraction B^Bh = 4 x 10~ 4 was assumed. 
The corresponding plot from the other approach is very similar. The h c signal is evident. 
The overall efficiencies determined from the Monte Carlo sample are 13.4% and 14.6% for 
the two inclusive approaches. Input values of M(h c ) and B^Bh are well reproduced. Results 
of Monte Carlo studies lead to the conclusion that the resonance fits to the data may be 
expected to have significance levels of ~4<r, statistical error on the mass of ~±0.6 MeV, and 
central values of the mass are reproduced within ~±0.6 MeV of the generated M(h c ). 

Figure 12(b) shows the data and the fit using the Monte Carlo optimized criteria for 
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TABLE II: Summary of estimated systematic uncertainties and their sums in quadrature. N/A 
means not applicable. 





M(h c 


) (MeV) 


B^B h 


x 10 4 


Systematic Uncertainty 


Inch 


Excl. 


Inch 


Excl. 


Number of i/j(2S) 


N/A 


N/A 


0.1 


N/A 


B(ib(2S) -» 7%) 


N/A 


N/A 


N/A 


0.8 


Background shape 


0.3 


0.2 


0.2 


0.3 


it® energy scale 


0.2 


~0 


0.1 


Signal shape 


0.1 


0.1 


0.3 


0.2 


h c width 


0.1 


0.1 


0.3 


0.2 


7T° efficiency 


~0 


~0 


0.1 


0.3 


photon efficiency 


~0 


~0 


0.2 


0.2 


Binning, fitting range 


0.1 


0.1 


0.3 


0.2 


Modeling of h c decays 


0.1 


0.3 


0.3 


~0 


r\ c mass 


0.1 


0.2 


0.1 


0.1 


r\ c width 


~0 


~0 


0.2 


0.1 


r] c branching ratios 


N/A 


~ 


N/A 


0.1 


Sum in quadrature 


±0.4 


±0.5 


±0.7 


±1.0 



the same inclusive approach as in Fig Efa)- Features in the Monte Carlo such as signal 
width, signal to background ratio, and approximate background shape mirror the data 
faithfully. The recoil spectrum and the fit for the other inclusive approach are very similar. 
Fit significance is approximately 3.8<r. Results from the two inclusive approaches differ by 
small amounts, with differences from the averages in M(h c ) of ±0.5 MeV and in B^Bh of 
±0.05 xlO -4 . The average results are listed in Table QJ 

The h c yield from the recoil mass against 7r° in the inclusive analysis is studied as a 
function of the angular distribution of the h c — > ^r\ c photon. The h c yield, shown in Fig.|3fa), 
is found to follow a l±cos 2 9 distribution (x 2 /d.o.f. = 1.7/2) as expected for an El transition 
from a spin 1 state. The background yield, shown in Fig. Efb), is uniform in cos#. The 
h c yield in the exclusive analysis is not sufficient to draw any conclusions regarding the 
corresponding angular distribution. 

Systematic uncertainties in the two analyses due to various possible sources have been 
estimated. Many sources are common, such as choice of background parameterization, h c 
resonance intrinsic width (r = 0.5 — 1.5 MeV), tt° line shape, bin size, and fitting range. The 
uncertainty in the branching ratio for ip(2S) —>■ ^rjc enters the systematic uncertainty for the 
exclusive analysis only while the uncertainty on the number oiijj(2S) decays applies to the 
inclusive analysis only. The estimated contributions are listed in Table ILT1 For the inclusive 
(exclusive) analysis they sum in quadrature to ±0.4 (0.5) MeV in M(h c ) and ±0.7 (1.0) x 10~ 4 
in B^Bh- The largest systematic error for the exclusive analysis, B(ip(2S) — > 7r/ c ), cancels 
in the ratio and we obtain B^B h /B v = 0.178 ± 0.049 ± 0.018. 

To summarize, we have observed the h c state, the l P\ state of charmonium, in the reaction 
ip(2S) — > 7T°h c , h c — > '-frjc, in exclusive and inclusive analyses. The significance of our 
observation is greater than 5a under a variety of methods to evaluate this quantity. We 



combine the results of the exclusive and inclusive analyses to obtain M(h c ) = 3524.4 ± 0.6 ± 
0.4 MeV and B(ip(2S) -> n°h c ) x B(h c -> jr) c ) = (4.0 ± 0.8 ± 0.7) x 10~ 4 . The following 
value is obtained for the hyperfine splitting: 

AM ?i/ ((M( 3 P J ))-M( 1 P 1 )) = +1.0 ± 0.6 ± 0.4 MeV. 

Thus, the combined result for M(h c ) is consistent with the spin-weighted average of the Xcj 
states. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science 
Foundation and the U.S. Department of Energy. 
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